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Abstract 

High energy protons produced by various sources of cosmic rays, e.g. supernova 
remnants, pulsar wind nebulae, active galactic nuclei and gamma-ray bursts participate 
in pj and pp interactions. Although pp interactions may be the dominant mechanism 
in our Galaxy, it is unclear how important p7 process is. We show that the upper 
bound on the fraction of protons participating in p"f interactions inside all Galactic 
astrophysical sources of cosmic rays is 10%. 
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r5 1 Introduction: 

High energy protons are expected to be emitted by various astrophysical objects e.g. su- 
\Q [ pernova remnants (SNRs), pulsar wind nebulae (PWN), microquasars, magnetars, active 

galactic nuclei (AGN) and gamma-ray bursts (GRBs). The relativistic protons interact with 
other protons and photons. High energy neutrinos and photons are subsequently produced 
in these interactions. Protons and other charged particles are deflected by Galactic and 
extra-galactic magnetic field before they arrive us. It is possible to trace the sources of 
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cosmic rays by detecting neutral particles e.g. neutrinos, photons from them. High energy 
neutrinos cannot be produced inside astrophysical sources of cosmic rays if they contain only 
leptons. Hence, detection of high energy neutrinos from a source would also be a signature of 
the existence of baryons inside that source. A large number of neutrino telescopes are under 
operation/construction to explore the neutrino universe in the near future. Innovative detec- 
tion methods have been employed in building these neutrino detectors PQEHHHHEIIEIIZIIH]- 
AMANDA, the high energy muon and neutrino detector in operation near the south pole has 
put limits on the neutrino fluxes from various point sources [HQ] in the northern hemisphere. 
The current limit on the diffuse all flavor neutrino flux is 3.3 x 10~ 7 GeVcm _2 sec _1 sr _1 in the 
energy range of 16TeV to 2PeV [ID] . This limit is from the data taken by AMANDA-II be- 
tween 2000 to 2003. In near future a more conservative limit of 4.2 x 10 _9 GeVcm _2 sec _1 sr _1 
or a real detection of the neutrino background is expected from IceCube experiment after 
three years of operation [IT] . 

Hadrons and leptons are expected to be shock accelerated inside astrophysical objects 
by the Fermi mechanism to very high energies. Charged and neutral pions can be produced 
in proton photon, proton proton or nucleon nucleon interactions, and the charged pions 
subsequently decay to produce neutrinos. Neutrino production interactions are py —> 7r + n — > 
e + v^v e v^n and pp — > 7r ± —>■ e^v^vJ^D^v^. We assume pn interactions are similar to pp 
interactions. The neutral pions decay to 7-rays. It is difficult to know what fractions of all 
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protons in our Galaxy participate in photo-pion (jyy) production and pp interactions. The 
positrons lose energy during their propagation unlike neutrinos. The propagation of electrons 
and positrons in our Galaxy has been studied earlier [12] . It has been suggested |13j that the 
detection of Glashow resonance [2] v e events at 6.3PeV by IceCube detector may offer us 
an opportunity to know whether pp interactions are the dominant mechanisms of ultrahigh 
energy neutrino production in the universe. The ratios of v e to muon flavor neutrinos (y e : 
{y^ + Ufj)) at the Glashow resonance energy (6.3 PeV) would be different for neutrinos 
produced in p'j and pp interactions, and by measuring these ratios it would also be possible 
to estimate neutrino mixing angles at ultrahigh energies |15j . However, effects like energy loss 
by muons is expected to reduce [IE] and the decay of pT and /i + produced in annihilation 
of secondary photons is expected to enhance the flux of v e [17]. These effects make the 
discrimination between pry and pp sources of neutrinos, by detecting Glashow resonance 
events, very difficult. Exploring other possibilities to reveal the underlying mechanisms of 
neutrino production in astrophysical sources is still of great interest. In our Galaxy the major 
sources of cosmic rays are SNRs and PWN. High energy gamma rays have been detected 
from many SNRs and PWN [18] although the leading mechanism (hadronic or leptonic) of 
7 ray production in these sources is not yet known. They are also expected to be sources 
of high energy neutrinos [19]. Protons and other heavy nuclei are accelerated in diffusive 
shocks. They interact with matter and radiation, lose energies and some of them escape from 
the sources, pp and pj interactions are the underlying mechanisms of neutrino production 
in these sources as discussed earlier. The gamma ray flux detected by CANGAROO [2E] 
from SN 1006 has been fitted with the theoretical model of it decay [21J. The flux produced 
in inverse Compton scattering of electrons by low energy photons is much lower than that 
produced in 7r° decay. Recently, the evidence of hadronic origin of gamma rays from Vela 
X [22] and PSR B 1259-63 [23] has been reported in two other papers. If a gamma ray 
energy flux of 2 x 10 -13 erg cm~ 2 s _1 is deteced from SN 1987A in the energy range of 0.1 
to lOTeV then it would imply that these gamma rays are of hadronic origin produced in tt° 
decay [21] . The neutrino fluxes expected from various PWN were calculated [25J considering 
detailed modeling of the sources and found to be detectable in some cases by large scale 
neutrino telescopes of km 2 area. High energy neutrino emission from py process inside the 
magnetospheres of isolated neutron stars have been discussed for magnetars [26] and for 
young pulsars [21]. In this paper the subject of our discussion is whether it is possible to 
constrain the parameter x, which denotes the fraction of the total number of Galactic protons 
participating in p7 interactions. 



2 Diffuse electron and positron spectra: 

Since high energy electrons and positrons lose energies faster than other charged particles, 
by various energy loss mechanisms e.g. bremsstrahlung, synchrotron emission and inverse 
Compton scattering, their observed fluxes are most likely of Galactic origin. Above 1 GeV the 
dominant energy loss processes are inverse Compton and synchrotron. It has been discussed 
in [2H] that approximately 10% of electrons and all positrons are secondary particles produced 
in nuclear interactions. There are some discrepancies in the measured data on electrons and 
positrons. Above several GeV, the large spread in the data can not be explained as an effect 
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of solar modulation. Systematic errors in the measurements might be responsible for the 
disagreement in the data taken by different experimental groups. Miiller [29J and Du Vernois 
et al. [30] suggested rescaling of the upper set of data points to match the electron flux at 
20 GeV obtained with the best fit of the lower set of data points that include the recent 
results of the AMS experiment [31] . After rescaling, all the data points in the energy range 
of 3 GeV to 2 TeV can be fitted with a single power law with the spectral index 3.44 ± 0.03 
[32] . With a similar procedure of renormalization of the data from MASS 89 [33] and MASS 
91 [34] experiments, the positron data points can be fitted with a single power law with the 
spectral index 3.43 ± 0.05 above 0.7 GeV. The single power law fitted spectra are shown in 
Fig. 4. and Fig. 5. of [32J, and we have used their results in the present paper. 

3 p7, pp interactions and e ± fluxes: 

Cosmic rays with energy up to 10 18 eV are most likely of Galactic origin. Protons can 
be accelerated to very high energies inside the Galactic cosmic ray sources, e.g. SNRs, 
PWN, and microquasars. One interesting feature is that interactions of shock accelerated 
protons with photons (jry — > e + Vpp ll v e n) leading to the production of neutrinos generate 
high energy positrons, but not electrons. Most other processes, on the other hand, produce 
equal number of electrons and positrons through pairs. These processes include py — * pe + e~ , 
pp — > e ± v i] y il v e {p ( ^) and 77 — > e + e~ . Such an asymmetry makes it possible to estimate the 
contribution of py process through evaluating the measured e + and e~ fluxes. 

The shock accelerated electron and proton spectra inside the sources could be approx- 
imated as power laws in energy with similar spectral index, i.e. S e -(E e ) = A e E~ a and 
Sp(Ep) = A p E~ a . The possible values of a can be calculated by numerical simulations 
[55] . Electrons and positrons lose energy faster than protons inside the sources. We have 
discussed in section 5 of this paper how the internal magnetic field and size of a source 
determine the escape probability of electrons and positrons from that source. We denote the 
electron and positron spectra injected into the Galaxy as Q e -(E e ) and Q e +(E e ) respectively. 
During propagation the cosmic ray protons and nucleons would mainly interact with Galactic 
matter and produce secondary electrons and positrons, since the Galactic diffuse gamma-ray 
background is not high enough for significant positron production by the propagating cosmic 
rays through pj interactions at the A resonance. The solution of the transport equation for 
the propagation of in the Galaxy is N(E e ) ~ Q(E e )V source T(E e ) /V occ (E e ) , where Q(E e ) 
is the source flux, V source is the volume of the sources, V occ (E e ) is the volume in the Galaxy 
occupied by electrons or positrons of energy E e , and r(E e ) is the lifetime of the particle with 
energy E e [2E] ■ We assume that out of the total number of protons inside all Galactic astro- 
physical objects, a fraction x participates in py interactions to produce neutrinos, another 
fraction y participates in pp interactions and the rest (z) escapes without any interactions 
or produce e ± pairs in py interactions, hence, x + y + z = 1. The injection source fluxes of 
electrons and positrons can be expressed as 



Q e -(E e )=P e -(E e ) + E(E e ) , 



(1) 



Q e +(E e ) = xf4> P7 ^ u+e +(E e ) + F(E e ) . 



(2) 
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Here F(E e ) denotes the secondary electron or positron flux injected into the Galaxy from 
all the interactions in which electron and positron pairs are produced, including inside the 
astrophysical sources and also during the propagation of the protons and nucleons in the 
Galaxy. Rotation-powered pulsars inject electron-positron pairs from the magnetosphere 
through magnetic one-photon (jB — > e + e~) pair production processes [36l ETJ [38] . WIMP 
annihilations also generate electron-positron pairs [39]. All these processes are included in 
our F(E e ) terms. We assume that astrophysical sources contain protons and electrons, but 
positrons are produced only in interactions. Thus we denote P e - (E e ) as the primary diffuse 
electron flux directly injected from all the Galactic astrophysical objects (not generated 
through interactions). In eqn(j2]), 4> P1 ^ v+e +(E e ) is the positron flux produced in photo- 
pion interactions normalized to one pion per nucleon (jry — > 7r + n — > e + v ll b> e iJ ll n) inside all 
astrophysical sources, assuming that all protons participate in pj interactions to produce 
neutrinos. The positrons produced in this way lose energies by bremsstrahlung, synchrotron 
and inverse Compton mechanisms before their escape from the sources. The emitted fluxes 
are related to the fluxes produced inside the sources by a factor / < 1. This factor accounts 
for the escaping probability and energy losses of electrons and positrons. The factor x denotes 
the fraction of protons which participate in pj interactions to produce neutrinos. If x = 
there is no neutrino production through pry interactions inside the sources, and x — 1 if all 
the neutrinos are produced in pry interactions. One can than relate the observed and source 
fluxes of electrons and positrons through the following ratios, 

N e -(E e ) -N e+ {E e ) = Q e -{E e )-Q e+ {E e ) _ P e -{E e ) - xf^ v+e+ {E e ) ^ 

N e -(E e ) Qe-(Ee) 

Besides pry interactions, all other interactions discussed earlier e.g. pp, 77, 'yB, WIMP 
annihilations, generate equal fluxes of electrons and positrons. The observed positron flux 
is about 10% of the observed electron flux and they are expected to be produced in various 
interactions (pp, py etc.) as discussed earlier. Hence, it is reasonable to assume that the 
flux of secondary electrons produced in interactions inside astrophysical sources and inside 
the Galaxy during the propagation of cosmic rays is much less compared to their original 
flux which is not produced in any interaction [28]. From this argument it follows that 
F(E e ) << P e -(E e ). In pry interactions charged and neutral pions are produced with almost 
equal probabilities. Unlike positrons, the nucleons produced in pry interaction escape from the 
sources before losing their energies significantly. The charged pion carries about 20% of the 
initial proton's energy and the final state leptons share the energy of 7r + equally. The energy 
E p of the proton and the energy E e of the positron produced in one p7 — > n + n — > e + u e i/^u tJ/ n 
interaction are related as E e = | x \E V . The total energies carried by all the positrons of 
energy E e produced in py interactions and all the protons of energy E p are related as [40j, 

E 2 e <P P ^ v+e+ (Ee) = \ x \ x \ElA v E~ a (4) 

where A p E~ a is the shock accelerated proton number flux per unit energy (at dE p ) inside 
the sources of cosmic rays, and the factor 1/2 takes into account the contribution of protons 
that produce 7r° (and hence, not contributing to e + production). The positron flux produced 
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in p7 interactions can be then expressed as 



{E e 



Sip ±-i p 



2-a 



40 El 

Eqn([n]) can then be expressed in terms of spectral index a, lepton energy E e 



(5) 
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P e -(E e ) + F(E e ) l + r (E 



fA e E~* 
(6) 

where we have used P e -(E e ) + F(E e ) = Q e -(E e ) from Eqn.(TjQ) and the ratio of observed 
positron and electron fluxes [32j |JT] is r(E e ) = N e +(E e )/N e -(E e ). Q e -(E e ) can be greater 
than fS e - (E e ) (the diffues electron flux emitted from all Galactic sources) at the most by a 
factor of [l+r(i? e )]. We have replaced Q e - (E e ) with [l+r(E e )]f S e -(E e ) in the denominator of 
the second term on the right hand side of Eqn.fJH]). The positron's energy E e and the proton's 
energy E p are related as E e = E p /20. Assuming that the shock accelerated electrons of mass 
m e and protons of mass m p have the same values of spectral indices the ratio of their shock 

accelerated fluxes at relativistic energies (E > m p c 2 ) is = C(^)^ ^ [28j , where ( is 
the ratio of the non relativistic electron and proton fluxes at any energy Ej. < m e c 2 (for 
non relativistic particles it is the kinetic energy). The numbers of non relativistic primary 
electrons and protons inside the sources are expected to be equal hence, £ ~ 1. We can then 
constrain x using the following expression, 



x 
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1 + r(E e 



2-a 



a-l 
2 



r(E e 



F(E e 



P e -(E e ) + F(E e 



(7) 



The ratio of observed positron and electron fluxes r(E e ) is available from observed data. 
However, we do not know how much is the ratio F(E e ) / [P e - (E e )+F(E e )} . Assuming F(E e ) = 
0, one can put an upper bound on x. 



x < 40 



l+r{E e 



r 1 x 2 - Q /m f 



q-l 
2 



r(E e ) . 



(8) 



Thus it is possible to constrain x using the observed data on positron and electron fluxes. 



4 Results: 

The current data on cosmic ray positron flux is up to about 60 GeV. In j4T] the data are 
given upto 30 GeV. Below 10 GeV solar modulation effect becomes important hence we 
have considered the energy range above this energy. As the ratio of the measured fluxes 
of positrons and electrons is small the errors involved in current measurements of this ratio 
affect the result. In future with more precise measurements the flux ratio could be determined 
more accurately. The average value of the observed flux ratios is 0.1. We calculate the upper 
bound on x for an energy range of lOGeV to 60GeV using the power law fit to the compilation 
of data from [32]. Our result is x < 0.1 for a = 2.1 — 2.4. The upper limit on x slightly 
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decreases with increasing a. In general the factor x depends on energy, but in the 10-60 
GeV range it is insensitive to the energy if we use the power law fit to the compiled data 
from [32] . In future it would be possible to constrain x at higher energies when the positron 
fluxes at higher energies are available. 

5 Discussions and Conclusion: 

In order to contribute to the observed positron flux at earth, positrons must not cool before 
escaping the source. This requires that the cooling time scale t coo i is longer than the escaping 
time scale t esc . For standard synchrotron cooling, the cooling time scale can be estimated as 

tcool = 7 9~R2"> V*) 

where m e is mass of a positron, c is the speed of light, 7 e is Lorentz factor of the positron, 
B is the internal magnetic field of the source, and a? is Thomson cross section. To estimate 
the escaping time scale, we adopt the random walk approximation. In a random magnetic 
field, the mean free path for each scattering may be estimated as the Larmor radius Rl- 
The number of scattering of a positron before escaping the source is ~ ("j^) 2 ' where R is the 
radius of the particle accelerating region inside a source. We therefore have 

R \iRl 



*- = (£)f- < 10 > 



2 

Since Rl = me ^ , it is interesting to observe that t coo \ and t esc have the same j e dependence. 
This suggests that for a particular source, either all positrons escape the source (t coo i > t esc ) 
or all positrons cool before escaping {t coo i < t esc ). With eqn. (fT0l and eqn.Q, the positron- 
escaping condition is 

fi_ 3 < 3.4R" 2/3 (11) 

where the convention Qi = Q/10 1 is adopted in cgs units. The typical values of B required 
to explain the observed radio and x-ray data from SN 1006 and SN 1987A are about 150/iG 
[2T] and lOmG [24] , respectively. The size of the positron emission region may vary from 
sources to sources. We can see that cosmic sources parameters are marginal to satisfy the 
constraint. Generally we expect some sources can satisfy eqn. ( TlTi) and contribute to the 
observed diffuse positron spectra. The factor / therefore denotes the fraction of all the 
jry positron generation sources that satisfy eqn. ([11]). The upper bound on the fraction of 
Galactic protons participating in jry interactions in our Galaxy is found to be 10%. Our 
current result is relevant to the positron energy range between 10 GeV to 60 GeV. In future 
it would be possible to extend our results with more observational data on positron fluxes 
at higher energies. 
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